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We give an introduction to the heavy-quark effective theory and the l/rriQ expan- 
sion, which provide the modern framework for a systematic, model-independent 
description of the properties and decays of hadrons containing a heavy quark. We 
discuss the applications of these concepts to spectroscopy and to the weak decays 
of B mesons. 



1 Introduction 

The weak decays of hadrons containing a heavy quark are employed for tests 
of the Standard Model and measurements of its parameters. They offer the 
most direct way to determine the weak mixing angles, to test the unitarity of 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix, and to explore the physics 
of CP violation. At the same time, hadronic weak decays also serve as a probe 
of that part of strong-interaction phenomenology which is least understood: 
the confinement of quarks and giuons inside hadrons. 

The structure of weak interactions in the Standard Model is rather simple. 
Flavour-changing decays are mediated by the coupling of the charged current 
to the M^-boson field. At low energies, the charged-current interaction gives 
rise to local four-fermion couplings, whose strength is governed by the Fermi 
constant 

Gf = / = 1.16639(2) GeV-2 . (1) 

According to the structure of the these interactions, the weak decays of hadrons 
can be divided into three classes: leptonic decays, in which the quarks of the 
decaying hadron annihilate each other and only leptons appear in the final 
state; semileptonic decays, in which both leptons and hadrons appear in the 
final state; and non-leptonic decays, in which the final state consists of hadrons 
only. Representative examples of these three types of decays are shown in 
Fig. I 

The simple quark-line graphs shown in this figure are a gross oversim- 
plification, however. In the real world, quarks are confined inside hadrons, 
bound by the exchange of soft giuons. The simplicity of the weak interactions 
is overshadowed by the complexity of the strong interactions. A complicated 
interplay between the weak and strong forces characterizes the phenomenology 
of hadronic weak decays. As an example, a more realistic picture of a non- 
leptonic decay is shown in Fig. 0. Clearly, the complexity of strong-interaction 
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Figure 1: Examples of leptonic (B — > r i't), semileptonic (B" — » D+e v^), and non- 
leptonic {B*^ — > D+tt^) decays of B mesons. 



effects increases with the number of quarks appearing in the final state. Bound- 
state effects in leptonic decays can be lumped into a single parameter (a "de- 
cay constant"), while those in semileptonic decays are described by invariant 
form factors, depending on the momentum transfer between the hadrons. 
Approximate symmetries of the strong interactions help to constrain the prop- 
erties of these form factors. For non-leptonic decays, on the other hand, we are 
still far from having a quantitative understanding of strong-interaction effects 
even in the simplest decay modes. 

Over the last decade, a lot of information on heavy-quark decays has been 
collected in experiments at e^e~ and hadron colliders. This has led to a rather 
detailed knowledge of the flavour sector of the Standard Model and many of 
the parameters associated with it. Thfice have been several gff^t discover- 
ies in this field, such as B^-B^ mixing.BEl, charmless B decays HlI, and rare 
decays induced by penguin operatorsQu. The experimental progress in heavy- 
flavour physics has been accompanied by a significant progress in theory, which 
was related to the discovery of heavy-quark symmetry and the development 
of the heavy-quark effective theory (HQET). The excitement about these de- 
velopments is caused by the fact that they allow (some) model-independent 
predictions in an area in which "progress" in theory often meant nothing more 
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Figure 2: More realistic representation of a non-leptonic decay. 

than the construction of a new model, which could be used to estimate some 
strong-interaction hadronic matrix elements. In these notes, we explain the 
physical picture behind heavy-quark symmetry and discuss the construction, 
as well as simple applications, of the heavy-quark expansion. Because of lack 
of time, we will have to focus on some particularly important aspects, empha- 
sizing the main ideas and concepts of the HQET. A more complete discussion 
of the applications of this formalism to_heavy-flavour phenomenology can be 
found in some recent review articlea^EnJ. The reader is also encouraged to 
consult the earlier review papers t3~tJ on the subject. 

Hadronic bound states of a heavy quark with light constituents (quarks, 
antiquarks and gluons) are characterized by a large separation of mass scales: 
the heavy-quark mass mg is much larger than the mass scale Aqcd associated 
with the light degrees of freedom. Equivalently, the Compton wave length 
of the heavy quark (Ag l/mg) is much smaller than the size of the hadron 
containing the heavy quark (i?had ~ 1/Aqcd)- Our goal will be to separate the 
physics associated with these two scales, in such a way that all dependence on 
the heavy-quark mass becomes explicit. The framework injidiich to perform 
this separation is the operator product expansion (OPE) EaM. The HQET 
provides us with a convenient technical tool to construct the OPE. Before 
we start to explore in detail the details of this effective theory, however, we 
should mention two important reasons why it is desirable to separate short- 
and long-distance physics in the first place: 

• A technical reason is that after the separation of short- and long-distance 
phenomena we can actually calculate a big portion of the relevant physics 
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(i.e. all short-distance effects) using perturbation theory and renormali- 
zation-group techniques. In particular, in this way we will be able to 
control all logarithmic dependence on the heavy-quark mass. 

• An important physical reason is that, after the short-distance physics has 
been separated, it may happen that the long-distance physics simplifies 
due to the realization of approximate symmetries, which imply non-trivial 
relations between observables. 

The second point is particularly exciting, since it allows us to make statements 
beyond the range of applicability of perturbation theory. Notice that here 
we are not talking about symmetries of the full QCD Lagrangian, such as its 
local gauge symmetry, but approximate symmetries realized in a particular 
kinematic situation. In particular, we will find that an approximate spin- 
fiavour symmetry is realized in systems in which a single heavy quark interacts 
with light degrees of freedom by the exchange of soft gluons. 

At this point it is instructive to recall a more familiar example of how ap- 
proximate symmetries relate the long-distance physics of several observables. 
The strong interactions of pions are severely constrained by the approximate 
chiral symmetry of QCD. In a certain kinematic regime, where the momenta 
of the pions are much less than 1 GeV (the scale of chiral-symmetry breaking) , 
the long-distance physics of scattering amplitudes is encoded in a few "reduced 
matrix elements", such as the pion decay constant. An effective low-energy 
theory called chiral perturbation theory provides a systematic expansion of 
scattering amplitudes in powers of the pion momenta, and thus helps to de- 
rive the relations between different scattering amplitudes imposed by chiral 
symmetry 113. We will find that a similar situation holds for the case of heavy 
quarks. Heavy-quark symmetry implies that, in the limit where mg ^ Aqcd, 
the long-distance physics of several observables is encoded in few hadronic pa- 
rameters, which can be defined in terms of operator matrix elements in the 
HQET. 

2 Heavy-Quark Symmetry 

2.1 The Physical Picture 

There are several reasons why the strong interactions of systems contain- 
ing heavy quarks are easier to understand than those of systems containing 
only light quarks. The first is asymptotic freedom, the fact that the effective 
coupling constant of QCD becomes weak in processes with laiigeptnomentum 
transfer, corresponding to interactions at short-distance scales At large 
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distances, on the other hand, the couphng becomes strong, leading to non- 
perturbative phenomena such as the confinement of quarks and gluons on—a 
length scale i?had ~ 1/Aqcd ~ 1 fm, which determines the size of hadronsc3. 
Roughly speaking, Aqcd ^ 0.2 GeV is the energy scale that separates the 
regions of large and small coupling constant. When the mass of a quark Q 
is much larger than this scale, it is called a heavy quark. The quarks of the 
Standard Model fall naturally into two classes: up, down and strange are light 
quarks, whereas charm, bottom and top are heavy quarks.^ For heavy quarks, 
the effective coupling constant asimq) is small, implying that on length scales 
comparable to the Compton wavelength Ag ~ ^/ttlq the strong interactions 
are perturbative and similar to the electromagnetic interactions. In fact, the 
quarkonium systems {QQ), whose size is of order Xq/asimQ) <C i?had, sue 
very much hydrogen-like. 

Systems composed of a heavy quark and light constituents are more com- 
plicated, however. The size of such systems is determined by i?had; and the 
typical momenta exchanged between the heavy and light constituents are of 
order Aqcd. The heavy quark is surrounded by a most complicated, strongly 
interacting cloud of light quarks, antiquarks, and gluons. In this case it is the 
fact that Xq <C i?hadj i-e. that the Compton wavelength of the heavy quark is 
much smaller than the size of the hadron, which leads to simplifications. To 
resolve the quantum numbers of the heavy quark would require a hard probe; 
the soft gluons exchanged between the heavy quark and the light constituents 
can only resolve distances much larger than Ag. Therefore, the light degrees 
of freedom are blind to the flavour (mass) and spin orientation of the heavy 
quark. They experience only its colour field, which extends over large dis- 
tances because of confinement. In the rest frame of the heavy quark, it is in 
fact only the electric colour field that is important; relativistic effects such as 
colour magnetism vanish as mg — s- oo. Since the heavy-quark spin partici- 
pates in interactions only through such relativistic effects, it decouples. That 
the heavy-quark mass becomes irrelevant can be seen as follows: As mg — > oo, 
the heavy quark and the hadron that contains it have the same velocity. In the 
rest frame of the hadron, the heavy quark is at rest, too. The wave function 
of the light constituents follows from a solution of the field equations of QCD 
subject to the boundary condition of a static triplet source of colour at the 
location of the heavy quark. This boundary condition is independent of mg, 
and so is the solution for the configuration of the light constituents. 

It follows that, in the limit mg — s- oo, hadronic systems which differ only 
in the flavour or spin quantum numbers of the heavy quark have the same 

"Ironically, the top quark is of no relevance to our discussion here, since it is too heavy to 
form hadronic bound states before it decays. 
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configuration of tlicir liglit degrees of freedom Although this observa- 

tion still does not allow us to calculate what this configuration is, it provides 
relations between the properties of such particles as the heavy mesons -B, I?, 
B* and D* , or the heavy baryons A;, and Ac (to the extent that corrections 
to the infinite quark-mass limit are small in these systems). These relations 
result from some approximate symmetries of the effective strong interactions 
of heavy quarks at low energies. The configuration of light degrees of freedom 
in a hadron containing a single heavy quark with velocity v does not change 
if this quark is replaced by another heavy quark with different flavour or spin, 
but with the same velocity. Both heavy quarks lead to the same static colour 
field. For Nh heavy-quark flavours, there is thus an S\J{2Nfi) spin-flavour sym- 
metry group, under which the effective strong interactions are invariant. These 
symmetries are in close correspondence to familiar properties of atoms: The 
flavour symmetry is analogous to the fact that different isotopes have the same 
chemistry, since to a good approximation the wave function of the electrons 
is independent of the mass of the nucleus. The electrons only see the total 
nuclear charge. The spin symmetry is analogous to the fact that the hyperfine 
levels in atoms are nearly degenerate. The nuclear spin decouples in the limit 
me/m^v 0. 

Heavy-quark symmetry is an approximate symmetry, and corrections arise 
since the quark masses are not infinite. In many respects, it is complementary 
to chiral symmetry, which arises in the opposite limit of small quark masses. 
However, whereas chiral symmetry is a symmetry of the QCD Lagrangian in 
the limit of vanishing quark masses, heavy-quark symmetry is not a symmetry 
of the Lagrangian (not even an approximate one) , but rather a symmetry of an 
effective theory, which is a good approximation of QCD in a certain kinematic 
region. It is realized only in systems in which a heavy quark interacts predom- 
inantly by the exchange of soft gluons. In such systems the heavy quark is 
almost on shell; its momentum fluctuates around the mass shell by an amount 
of order Aqcd- The corresponding fluctuations in the velocity of the heavy 
quark vanish as Aqcd/'tiq 0. The velocitvJpecomes a conserved quantity 
and is no longer a dynamical degree of freedomEZl. Nevertheless, results derived 
on the basis of heavy-quark symmetry are model-independent consequences of 
QCD in a well-defined limit. The symmetry-breaking corrections can, at least 
in principle, be studied in a systematic way. A convenient framework for ana- 
lyzing these corrections is provided by the heavy-quark effective theory. Before 
presenting a detailed discussion of the formalism, we shall first point out some 
of the important implications of heavy-quark symmetry for the spectroscopy 
and weak decays of heavy hadrons. 
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2.2 Spectroscopic Implications 

The spin-flavour symmetry leads to many interesting relations between the 
properties of hadrons containing a heavy miark. The most direct consequences 
concern the spectroscopy of such states Ba. In the limit uiq — > oo, the spin 
of the heavy quark and the total angular momentum j of the light degrees of 
freedom inside a hadron are separately conserved by the strong interactions. 
Because of heavy-quark symmetry, the dynamics is independent of the spin and 
mass of the heavy quark. Hadronic states can thus be classified by the ouantum 
numbers (flavour, spin, parity, etc.) of the light degrees of freedom c3. The 
spin symmetry predicts that, for fixed j 7^ 0, there is a doublet of degenerate 
states with total spin J — j ± ^. The flavour symmetry relates the properties 
of states with different heavy-quark flavour. 

In general, the mass of a hadron Hq containing a heavy quark Q obeys an 
expansion of the form 

mH = mQ+A+^— + 0{l/ml). (2) 

The parameter A represents contributions arising froni_all terms in the La- 
grangian that are independent of the heavy-quark massE3, whereas the quan- 
tity Am^ originates from the terms of order 1 /mq in the effective Lagrangian 
of the HQET. For the moment, the detailed structure of these terms is of no 
relevance; it will be discussed at length in the next section. For the ground- 
state pseudoscalar and vector mesons, one can parametrize the contributions 
from thejymg corrections in terms of two quantities, Ai and A2, in such a 
way thatEil 

Am2 = -Ai + 2 [ J( J + 1) - 1] A2 . (3) 

Here J is the total spin of the meson. The first term, — Ai/2rnQ, arises from the 
kinetic energy of the heavy quark inside the meson; the second term describes 
the interaction of the heavy-quark spin with the gluon field. The hadronic pa- 
rameters A, Ai and A2 are independent of toq. They characterize the properties 
of the light constituents. 

Consider, as a first example, the SU(3) mass splittings for heavy mesons. 
The heavy-quark expansion predicts that 

TTiBs - "^Sd = As - Ad + 0(1 /mb) , 

rriDs - rriDa = A,, - A^ + 0(1/toc) , (4) 

where we have indicated that the value of the parameter A depends on the 
flavour of the light quark. Thus, to the extent that the charm and bottom 
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quarks can both be considered sufficiently heavy, the mass sphttings should 
be similar in the two systems. This prediction is confirmed experimentally, 
since t3 

rriBs - rriBa = (90 ± 3) MeV , 

TODs -mD^ = (99±1) MeV. (5) 

As a second example, consider the spin splittings between the ground-state 
pseudoscalar (J = 0) and vector (J = 1) mesons, which are the members of 
the spin-doublet with j = ^. The theory predicts that 

m%, - m% = 4A2 + 0{l/mb) , 

to|,. - m|, = 4A2 + 0(l/mc) . (6) 
The data are compatible with this: 

771%. - m| ~ 0.49 GeV^ , 

771%, - m|, ~ 0.55 GeV^ . (7) 

Assuming that the B system is close to the heavy-quark limit, we obtain the 
value 

A2 ~ 0.12 GeV^ (8) 

for one of the hadronic parameters in (|^). This quantity playSj^n important 
role in the phenomenology of inclusive decays of heavy hadronstl. 

A third example is provided by the mass splittings between the ground- 
state mesons and baryons containing a heavy quark. The HQET predicts that 

"^A6 - ms = Abaryon - A,ncson + 0{l/7nb) , 

rriA, - 777D = Abaryon " A,noson + 0(1/TOc) ■ (9) 

This is again consistent with the experimental results 
WA, - ms = (346 ± 6) MeV , 

TOA^ - 771D = (416 ± 1) MeV , (10) 

although in this case the data indicate sizeable symmetry-breaking corrections. 
For the mass of the A^ baryon, we have used the value 

mA, = (5625 ±6) MeV, (11) 

which is obtained by averaging the result El mA^ = (5639 ± 15) MeV with 
the value ttia, = (5623 ± 5 ± 4) MeV reported by the CDF CoUaborationy. 
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The dominant correction to the relations comes from the contribution of 
the chromo-magnetic interaction to the masses of the heavy mesons,^ which 
adds a term 3X2 /2mQ on the right-hand side. Including this term, we obtain 
the refined prediction that the values of the following two quantities should be 
close to each other: 

mA, -niB- = (312 ± 6) MeV , 

TOA^ -mi) - — ^ = (320± 1) MeV (12) 

This is clearly satisfied by the data. 

The mass formula can also be used to derive information on the heavy- 
quark (pole) masses from the observed hadron masses. Introducing the "spin- 
averaged" meson masses TOB = j (mB+imB*) — 5.31 GeVandmD = j (mu + 
3m£i.) ~ 1.97 GeV, we find that 



mb - TUc = {niB - mo) 1 1 - „— ^— l-0(l/mo)l, (13) 

1^ 2m,BmD ) 

where ©(l/mg) is used as a generic notation representing terms suppressed 
by three powers of the b- or c-quark maasea_. Using theoretical estimates for 
the parameter Ai, which lie in the rangeo t3 

Ai = -(0.3±0.2) GeV^ (14) 

this relation leads to 



TOfa-rric = (3.39 ±0.03 ±0.03) GeV, (15) 

where the first error reflects the uncertainty in the value of Ai, and the second 
one takes into account unknown higher-order corrections. 



2.3 Exclusive Semileptonic Decays 

Semileptonic decays of B mesons have received a lot of attention in recent 
years. The decay channel B D*£D has the largest branching fraction of all 
B-meson decay modes. From a theoretical point of view, semileptonic decays 
are simple enough to allow for a reliable, quantitative description. The analysis 
of these decays provides much information about the strong forces that bind 
the quarks and gluons into hadrons. Heavy-quark symmetry implies relations 

''Because of the spin symmetry, there is no such contribution to the masses of the Aq baryons. 
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between the weak decay form factors of heavy mesons, which ate of particular 
interest. These relations have been derived by Isgur and WiseE3, generalizing 
ideas developed by Nussinov and WetzelE3, and by Voloshin and ShifmanEJS. 

Consider the elastic scattering of a B meson, B{v) — > B{v'), induced by a 
vector current coupled to the b quark. Before the action of the current, the light 
degrees of freedom inside the B meson orbit around the heavy quark, which 
acts as a static source of colour. On average, the b quark and the B meson have 
the same velocity v. The action of the current is to replace instantaneously (at 
t = to) the colour source by one moving at a velocity v' , as indicated in Fig. |[ 
li V ~ v', nothing happens; the light degrees of freedom do not realize that 
there was a current acting on the heavy quark. If the velocities are different, 
however, the light constituents suddenly find themselves interacting with a 
moving colour source. Soft gluons have to be exchanged to rearrange them so 
as to form a B meson moving at velocity v' . This rearrangement leads to a 
form-factor suppression, which reflects the fact that as the velocities become 
more and more different, the probability for an elastic transition decreases. 
The important observation is that, in the limit mi, — > oo, the form factor can 
only depend on the Lorentz boost j — v-v' that connects the rest frames of the 
initial- and final-state mesons. Thus, in this limit a dimensionless probability 
function ^ (?;•«') describes the transition. It is called the Isgur-Wise functioncj. 
In the HQET, which provides the appropriate framework for taking the limit 
nih oo, the hadronic matrix element describing the scattering process can 
thus be written as 

— {B{v')\K,j''b, \B{v)) = ^{v ■ v') {v + vY . (16) 
niB 

Here, by and by' are the velocity-dependent heavy-quark fields of the HQET, 
whose precise definition will be discussed in Sec. ^. It is important that the 
function ^{v ■ v') does not depend on rrifc. The factor I/tob on the left-hand 
side compensates for a trivial dependence on the heavy-meson mass caused by 
the relativistic normalization of meson states, which is conventionally taken to 
be 

{B{p')\B{p)) ^2mBv'' {2nf 5\p- p') . (17) 

Note that there is no term proportional to {v — v')^ in (^6|). This can be seen 
by contracting the matrix element with (v — w')^, which must give zero since 
fby ~ by and bv'f = by'. 

It is more conventional to write the above matrix element in terms of an 
elastic form factor Fe\{q'^) depending on the momentum transfer q'^ — [p — p'Y- 

{B{v')\brb\B{v))^F,,{q^){p + pr . (18) 
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t<tQ t = tQ t>tQ 

Figure 3: Elastic transition induced by an external heavy-quark current. 

where p^'^ = msv^'^ . Comparing this with (p^, we find that 

F,,{q') ^ i{vv') , ^-2ml{vv' ^1). (19) 

Because of current conservation, the elastic form factor is normahzed to unity 
at = 0. This condition imphes the normalization of the Isgur-Wise function 
at the kinematic point v ■ v' ^ \, i.e. for v = v': 

e(l) = l. (20) 

It is in accordance with the intuitive argument that the probabihty for an 
elastic transition is unity if there is no velocity change. Since for v — v' the 
daughter meson is at rest in the rest frame of the parent meson, the point 
V ■ v' = \ is referred to as the zero-recoil hmit. 

We can now use the flavour symmetry to replace the b quark in the final- 
state meson by a c quark, thereby turning the B meson into a D meson. Then 
the scattering process turns into a weak decay process. In the infinite mass 
limit, the replacement h^i — > c^i is a symmetry transformation, under which 
the effective Lagrangian is invariant. Hence, the matrix element 

, ^ {D{v')\ c,,j%, \B{v)) - ■ v') {v + v'r (21) 

is still determined by the same function ^{v ■ v'). This is interesting, since in 
general the matrix element of a flavour-changing current between two pseu- 
doscalar mesons is described by two form factors: 

{D{v')\ -c^h \B{v)) = U{q') {p + p'Y - /-(?') (P - P'Y ■ (22) 

Comparing the above two equations, we flnd that 

2n mB± rriD 
2^mBmD 

q^ = + m?jj — 2mBniD v ■ v' . (23) 
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Thus, the heavy-quark flavour symmetry relates two a priori independent form 
factors to one and the same function. Moreover, the normahzation of the 
Isgur-Wise function at v ■ v' — 1 now imphes a non-trivial normalization of 



the form factors f±{q^) 



at the point of maximum momentum transfer, 



/±(9max) 



(24) 



The heavy-quark spin symmetry leads to additional relations among weak 
decay form factors. It can be used to relate matrix elements involving vector 
mesons to those involving pseudoscalar mesons. A vector meson with longi- 
tudinal polarization is related to a pseudoscalar meson by a rotation of the 



heavy-quark spin. Hence, the spin-symmetry transformation c|,, 



B ^ D with B D* transitions. The result of this transformation is El: 



', jslates 



1 



1 



{D*iv',e)\c,,rb,\Biv)) = ie'^-"'' £»0 ■ v') , 



^riiBniD* 



{D*{v\e)\c„,Y'l5bv\B{v)) = 6*^" {vv' + 1) ~ v'^' e* ■ v ^{v ■ v') 



(25) 



where e denotes the polarization vector of the D* meson. Once again, the 
matrix elements are completely described in terms of the Isgur-Wise function. 
Now this is even more remarkable, since in general four form factors, V{q'^) 
for the vector current, and Ai{q'^), i = 0, 1, 2, for the axial vector current, are 
required to parametrize these matrix elements. In the heavy-quark limit, they 
obey the relations EZl 



± Tod. 
2y/mBmD* 



1 - 



{rriB + muY 



Ai{q^). 



q^ = m\ + m^* — 2mBiTiD- v 



(26) 



Equations ( P3| ) and ( p6| ) summarize the relations imposed by heavy-quark 
symmetry on the weak decay form factors describing the semileptonic decay 
processes B ^ Dlv and B D*tv. These relations are model-independent 
consequences of QCD in the limit where m;,, TOc 3> Aqcd- They play a crucial 
role in the determination of the CKM matrix element | Vcb \ ■ In terms of the 
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recoil variable w — V-^v' , the differential semileptonic decay rates in the heavy- 
quark limit become t 



dT{B Dlv) 
dw 

dr(B D*iv) 



G p 
487r3 

Jjj iKbP (ms - mc. ml, ^w^l {w + if 



4w m'i, — 2w mBmu* + m^j 



w + 1 



{ms ~ mo'Y 



These expressions receive symmetry-breaking corrections, since the masses of 
the heavy quarks are not infinitely heavy. Perturbative corrections of order 
a"(mQ) can be calculated order by order in perturbation theory. A more 
difficult task is to control the non-perturbative power corrections of order 
( Aqcd /mq ) " . The HQET provides a systematic framework for analysing these 
corrections. For the case of weak-decay form factors, the analysis of the l/mg 
corrections was performed by LukeEj. Later, Falk and the present author have 
also analysed the structure of corrections for both meson and baryon 

weak decay form factors 0. We shall not discuss these rather technical issues 
in detail, but only mention the most important result of Luke's analysis. It 
concerns the zero- recoil limit, where an analoeup of the Ademollo-Gatto thco- 
remllJ can be proved. This is Luke's thcoremEj, which states that the matrix 
elements describing the leading I/toq corrections to weak decay amplitudes 
vanish_aLjzero recoil. This theorem is valid to all orders in perturbation the- 
ory tiloH Most importantly, it protects the B —±^*£D decay rate from 
receiving first-order l/mq corrections at zero recoil (A similar statement 
is not true for the decay B ^ D £ D, however. The reason is simple but some- 
what subtle. Luke's theorem protects only those form factors not multiplied by 
kinematic factors that vanish for v — v' . By angular momentum conservation, 
the two pseudoscalar mesons in the decay B D £p must be in a relative p 
wave, and hence the amplitude is proportional to the velocity \vd\ of the D 
meson in the _B-meson rest frame. This leads to a factor (w^ — 1) in the decay 
rate. In such a situation, form factors that are kinematically suppressed can 
contributeEZl.) 



2.4 Model-Independent Determination of \Vcb\ 

We will now discuss the most important application of the HQET in the context 
of semileptonic decays of B mesons. A model-independent determination of the 
CKM matrix element \ Vcb \ based on heavy-quark symmetry can be obtained by 
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measuring the recoil spectrum of D* mesons produced in i? ^ D*£ v decays m. 
In the heavy-quark Umit, the differential decay rate for this process has been 
given in (^7|). In order to allow for corrections to that limit, we write 



Aw 



(niB - mo* f m%, \J - \ {w + 1)^ 



1 + 



Aw mi, — 2w mBiriD* + mj 



1 



\Vcb\'T\w), 
(28) 



where the hadronic form factor T{w) coincides with the Isgur-Wise function 
up to symmetry-breaking corrections of order as{mQ) and Aqcd/^t-q- The 
idea is to measure the product \Vcb\ ^{w) as a function of w, and to extract 
\Vcb\ from an extrapolation of the data to the zero- recoil point w = 1, where 
the B and the D* mesons have a common rest frame. At this kinematic point, 
heavy-quark symmetry helps to calculate the normalization !F{1) with small 
and controlled theoretical errors. Since the range of w values accessible in this 
decay is rather small (1 < w < 1.5), the extrapolation can be done using an 
expansion around w — 1: 



T{w) = l-f{w-l) + 



(29) 



The slope is treated as a fit parameter. 

Measurements of JJie recoil spectrum have been performed first by the 
ARGUS El and CLEOll^ Collaborations in experiments operatiiig at the T(4s) 
resonance, and more recently by the ALEPHcJ and DELPHIcS Collaborations 
at LEP. As an example, Fig. ^ shows the data reported by the CLEG Collabo- 
ration. The results obtained by the various experimental groups from a linear 
fit to their data are summarized in Table |l|. The weighted average of these 
results is 



\Vcb\T{l) = (34.6 ±1.7) X 10" 
= 0.82 ± 0.09 . 



(30) 



The effiect of a positive curvature of the form factor has been investigated by 
Stone cIl, who finds that the value of \Vcb \ ^{^) may change by up to -t-4%. We 
thus increase the above value by (2 ± 2)% and quote the final result as 



\Vcb\J'{l) = (35.3 ± 1.8) X 10 



-3 



(31) 
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Figure 4: CLEO jifita for the product \Vcb\^iw), ^ extracted from the recoil spectrum in 
B — ► D*£u decaysEj. The hue shows a linear fit to the data. 



In future analyses, the extrapolation to zero recoil should be performed in- 
cluding higher-order terms in the expansion (^9|) . It can be shown in a model- 
independent way that the shape of the—facrn factor is highly constrained by 
analyticity and unitarity requirements c3'e3. In particular, the curvature at 
ui = 1 is strongly correlated with the slope of the form factor. For the value of 
'(P given in (|30|), one obtains a small positive curvature cj, in agreement with 
the assumption made in Ref. 47. 



Table f: Values for |V'c6|.^(l) (in units of 10 and ^ extracted from measurements of 
the recoil spectrum in B — > D*tv decays 





1^61^(1) (10-3) 




ARGUS 
CLEO 
ALEPH 
DELPHI 


38.8 ±4.3 ±2.5 
35.1 ±1.9 ±2.0 
31.4 ±2.3 ±2.5 
35.0 ±1.9 ±2.3 


1.17 ±0.22 ±0.06 
0.84 ±0.12 ±0.08 
0.39 ±0.21 ±0.12 
0.81 ±0.16 ±0.10 



Heavy-quark symmetry implies that the general structure of the symmetry- 
breaking corrections to the form factor at zero recoil isES 



= 1 + X 



A 



QCD 



const X 



^^QCD 



?7a(1 + <5i/™2), (32) 
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where tja is a short-distance correction arising from the (finite) renormahzation 
of the flavour-changing axial current at zero recoil, and (Si/m2 parametrizes 
second-order (and higher) power corrections. The absence of£rst-order power 
corrections at zero recoil is a consequence of LukelalheoremLJ. The one-loop 
expression for tja has been known for a long timeE3Ej'E2l: 

= 1 + ^ ( H^^LlHh In^^-l)^ 0.96. (33) 

The scale M in the running coupling constant can be_fixed by adopting the pre- 
scription of Brodsky, Lepage and Mackenzie (BLM)lj, according to which it is 
identified with the average virtuality of the gluon in the one-loop diagrams that 
contribute to rjA- If o^iM) is defined in the modified minimal subtraction (ms) 
scheme, the result isE3 M ~ O.bl^/mimb. Several estimates of higher-order 
corrections to r]A have been discussed. The next-to-leadipg order resumma- 
tion of logarithms of the type [ug \n.{mi, / nic)]" leads toE3o tja ~ 0.985. On 
the other hand, the resummation of "renormalon-chain" contributions of the 
form /3g~^Q;", where /3o is the first coefficient of the QCD /3- function, gives Ea 
r]A — 0.945. Using these partial resummations to estimate the uncertainty 
results in tja = 0.965 ± 0.020. Recently, Czarnecki has improved this estimate 
by calculating tja at two-loop orderES. His result, 

r;^ = 0.960 ± 0.007 , (34) 

is in excellent agreement with the ELM- improved one- loop estimate (|3^). Here 
the error is taken to be the size of the two-loop correction. 

The analysis of the power corrections i5i/m2 is more difficult, since it cannot 
rely on perturbation theory. Three approaches have been discussed: in the 
"exclusive approach" , all 1/ TOq operatprsrn the HQET are classified and their 
matrix elements estimated, leading toEMZl Si/m^ = —(3 ± 2)%; the "inclusive 
approach" has been used to derive the bound Si/m^ < —3%, and to estimate 
that E3'| (5i/™2 = -(7 ±3)%; the "hybrid approach" combines the virtues of 
the_former two to obtain a more restrictive lower bound on 6i/,j^2. This leads 
toB 

f^i/m^ = -0.055 ± 0.025. (35) 

Combining the above results, adding the theoretical errors linearly to be 
conservative, gives 

J^(l) = 0.91 ±0.03 (36) 

for the normalization of the hadronic form factor at zero recoil. Thus, the 
corrections to the heavy-quark limit amount to a moderate decrease of the 

'^This bound has been criticised in Ref. 59. 



16 



form factor of about 10%. This can be used to extract from the experimental 
result (^ij) the model-independent value 

iKbl = (38.8 ± 2.0exp ± 1.2th) X 10-3 . (37) 

3 Heavy-Quark Effective Theory 

3.1 The Effective Lagrangian 

The effects of a very heavy particle often become irrelevant at low energies. It 
is then useful to construct a low-energy effective theory, in which this heavy 
particle no longer appears. Eventually, this effective theory will be easier 
to deal with than the full theory. A familiar example is Fermi's theory of 
the weak interactions. For the description of weak decays of hadrons, the 
weak interactions can be approximated by point-like four-fermion couplings, 
governed by a dimensionful coupling constant Gp. Only at energies much 
larger than the masses of hadrons can the effects of the intermediate vector 
bosons, W and Z, be resolved. 

The process of pemoving the degrees of freedom of a heavy particle involves 
the following steps E^H^: one first identifies the heavy-particle fields and "in- 
tegrates them out" in the generating functional of the Green functions of the 
theory. This is possible since at low energies the heavy particle does not appear 
as an external state. However, although the action of the full theory is usually 
a local one, what results after this first step is a non-local effective action. The 
non-locality is related to the fact that in the full theory the heavy particle 
with mass M can appear in virtual processes and propagate over a short but 
finite distance Aa; ~ 1/M. Thus, a second step is required to obtain a local 
effective Lagrangian: the non-local effective action is rewritteii_aii_an infinite 
series of local terms in an Operator Product Expansion (OPE) Roughly 
speaking, this corresponds to an expansion in powers of 1/M. It is in this step 
that the short- and long-distance physics is disentangled. The long-distance 
physics corresponds to interactions at low energies and is the same in the full 
and the effective theory. But short-distance effects arising from quantum cor- 
rections involving large virtual momenta (of order M) are not reproduced in 
the effective theory, once the heavy particle has been integrated out. In a third 
step, they have to be added in a perturbative way using renormalization-group 
techniques. These short-distance effects lead to a renormalization of the coef- 
ficients of the local operators in the effective Lagrangian. An example is the 
effective Lagrangian for non-leptonic weak decays, in which radiative correc- 
tions from hard gluons with virtual momenta in the range between mw and 
some renormalization scale ~ 1 GeV give rise to Wilson coefficients, which 
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Figure 5: Philosophy of the heavy-quark effective theory. 



The heavy-quark effective theory (HQET) is constructed to provide a sim- 
plified description of processes where a heavy quark interacts-with light degrees 
of freedom predominantly by the exchange of soft gluons lZI"lZI. Clearly, mq 
is the high-energy scale in this case, and Aqcd is the scale of the hadronic 
physics we are interested in. The situation is illustrated in Fig. ^. At short 
distances, i.e. for energy scales larger than the heavy-quark mass, the physics is 
perturbative and described by ordinary QCD. For mass scales much below the 
heavy-quark mass, the physics is complicated and non-perturbative because 
of confinement. Our goal is to obtain a simplified description in this region 
using an effective field theory. To separate short- and long-distance effects, we 
introduce a separation scale /i such that Aqcd ^ M <C mq . The HQET will be 
constructed in such a way that it is identical to QCD in the long-distance re- 
gion, i.e. for scales below /i. In the short-distance region, the effective theory is 
incomplete, however, since some high-momentum modes have been integrated 
out from the full theory. The fact that the physics must be independent of the 
arbitrary scale ji allows us to derive renormalization-group equations, which 
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we shall employ to deal with the short-distance effects in an efficient way. 

Compared with most effective theories, in which the degrees of freedom 
of a heavy particle are removed completely from the low-energy theory, the 
HQET is special in that its purpose is to describe the properties and decays of 
hadrons which do contain a heavy quark. Hence, it is not possible to remove 
the heavy quark completely from the effective theory. What is possible is to 
integrate out the "small components" in the full heavy-quark spinor, which 
describe the fluctuations around the mass shell. 

The starting point in the construction of the low-energy effective theory 
is the observation that a very heavy quark bound inside a hadron moves more 
or less with the hadron's velocity v, and is almost on shell. Its momentum can 
be written as 

p^Q = rnQvi" + , (38) 

where the components of the so-called residual momentum k are much smaller 

than toq. Note that w is a four- velocity, so that = 1. Interactions of the 
heavy quark with light degrees of freedom change the residual momentum by 
an amount of order Afc Aqcd, but the corresponding changes in the heavy- 
quark velocity vanish as Aqcd /i^-q ^ 0. In this situation, it is appropriate to 
introduce large- and small-component fields, hy and H-^, by 

h,{x) = e^™*"-^ P+ Q{x) , H,{x) = e^™«"-" P_ Q{x) , (39) 

where P+ and P_ are projection operators defined as 

P. = ^. (40) 

It follows that 

Q{x) = e-^Q"- [K (x) + Hy (x)] . (41) 

Because of the projection operators, the new fields satisfy fhy = hy and i/} Hy = 
—Hy. In the rest frame, i.e. for u'* = (1,0,0,0), hy corresponds to the upper 
two components of Q, while Hy corresponds to the lower ones. Whereas hy 
annihilates a heavy quark with velocity v, Hy creates a heavy antiquark with 
velocity v. 

In terms of the new fields, the QCD Lagrangian for a heavy quark takes 
the form 

Cq = Q{i]f) - mq) Q 

= hy iv- D hy — Hy (iv- D + 2mQ) Hy 

+ hyip±Hy+Hyip±hy, (42) 
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where D'^ = 13'' — v^v-D\s orthogonal to the heavy-quark velocity: v ■ D± = 0. 
In the rest frame, Dj = (0, D ) contains the spatial components of the covariant 
derivative. From (|2|), it is apparent that hy describes massless degrees of 
freedom, whereas Hy corresponds to fluctuations with twice the heavy-quark 
mass. These are the heavy degrees of freedom that will be eliminated in the 
construction of the effective theory. The fields are mixed by the presence of the 
third and fourth terms, which describe pair creation or annihilation of heavy 
quarks and antiquarks. As shown in the first diagram in Fig. ^ in a virtual 
process a heavy quark propagating forward in time can turn into an antiquark 
propagating backward in time, and then turn back into a quark. The energy 
of the intermediate quantum state hhH is larger than the energy of the initial 
heavy quark by at least 2niQ. Because of this large energy gap, the virtual 
quantum fluctuation can only propagate over a short distance Ax l/mq. 
On hadronic scales set by -Rhad = V^QCD, the process essentially looks like a 
local interaction of the form 

Kip^TT—iP^h^^ (43) 
2mQ 

where we have simply replaced the propagator for Hy by l/2mQ. A more 
correct treatment is to integrate out the small-component field Hy, thereby 
deriving a non-local effective action for the large-component field hy, which 
can then be expanded in terms of local operators. Before doing this, let us 
mention a second type of virtual corrections involving pair creation, namely 
heavy-quark loops. An example is shown in the second diagram in Fig. |[ 
Heavy-quark loops cannot be described in terms of the effective fields hy and 
Hy, since the quark velocities inside a loop are not conserved and are in no 
way related to hadron velocities. However, such short-distance processes are 
proportional to the small coupling constant as{mQ) and can be calculated in 
perturbation theory. They lead to corrections that are added onto the low- 
energy effective theory in the renormalization procedure to be discussed later. 




Figure 6: Virtual fluctuations involving pair creation of heavy quarks. In the first diagram, 
time flows to the right. 

On a classical level, the heavy degrees of freedom represented by Hy can 
be eliminated using the equation of motion. Taking the variation of the La- 
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grangian with respect to the field H^, we obtain 

{ivD + 2mQ)H, = ip^h,. (44) 
This equation can formally be solved to give 

Hv^^ ^-—f,^P±K, (45) 

showing that the small-component field is indeed of order l/mg. We can 
now insert this solution into (^2|) to obtain the "non-local effective Lagrangian" 

CcS = hyiv-D hy + hvi p± — -ip±hy. (46) 

2niQ + iv-U 

Clearly, the second term corresponds to the first class of virtual processes 
shown in Fig. ||. 

It is possible to derive this Lagrangian in a more elegant way by manipu- 
lating the genfirating functional for QCD Green's functions containing heavy- 
quark fields im. To this end, one starts from the field redefinition ( [4l| ) and 
couples the large-component fields to external sources p^. Green's functions 
with an arbitrary number of hy fields can be constructed by taking derivatives 
with respect to py. No sources are needed for the heavy degrees of freedom 
represented by Hy. The functional integral over these fields is Gaussian and 
can be performed explicitly, leading to the effective action 



S'cff = y d^x/Iefr-ilnA, (47) 

with CcS as given in (^). The appearance of the logarithm of the determinant 

= eyi-p (^Tt \n[2mQ + iv D - iri\^ (48) 

is a quantum effect not present in the classical derivation presented above. 
However, in this case the determinant can be regulated in a gauge-invariant 
way, and by choosinELlhe axial gauge v ■ A = Q one shows that In A is just an 
irrelevant constant EZliIa. 

Because of the phase factor in (p]), the x dependence of the effective 
heavy-quark field hy is weak. In momentum space, derivatives acting on hy 
correspond to powers of the residual momentum k, which by construction is 
much smaller than toq. Hence, the non-local effective Lagrangian ( p6| ) allows 
for a derivative expansion in powers of iD/mQ: 

CcS = hyivDhy + ^ hyip^i ~'^-^^-—] ip^hy. (49) 
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Taking into account that /i„ contains a P+ projection operator, and using the 
identity 

2 



(50) 



where [iD^^ , iD'^] — igsG^" is the gluon field-strength tensor, one finds that 
Cff = K iv-D K + TT— {tD± f K + K fT^. G^" h, + 0{l/m^n) • (51) 



In the limit mg oo, only the first terms remains: 



(52) 



This is the effective Lagrangian of the HQET. It gives rise to the Feynman 
rules depicted in Fig. ^ 
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Figure 7: Feynman rules of the HQET and a are colour indices). A heavy quark is 
represented by a double line labelled by the velocity v and the residual momentum k. The 
velocity v is conserved by the strong interactions. 

Let us take a moment to study the symmetries of this Lagrangian^Zl Since 
there appear no Dirac matrices, interactions of the heavy quark with gluons 
leave its spin unchanged. Associated with this is an SU(2) symmetry group, 
under which C^o is invariant. The action of this symmetry on the heavy-quark 
fields becomes most transparent in the rest frame, where the generators S"* of 
SU(2) can be chosen as 



2 I 







(53) 



Here cr* are the Pauli matrices. An infinitesimal SU(2) transformation 
(1 + ie • S" ) leaves the Lagrangian invariant: 



5 Loo = hv [iv-D, ie-S]hv = . 



(54) 
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Another symmetry of the HQET arises since the mass of the heavy quark does 
not appear in the effective Lagrangian. For heavy quarks moving at the 
same velocity, eq. (|5^ ) can be extended by writing 

1=1 

This is invariant under rotations in flavour space. When combined with the 
spin symmetry, the symmetry group-jis, promoted to SU(2iV?i). This is the 
heavy-quark spin-flavour symmetry EaO. Its physical content is that, in the 
limit ruQ 00, the strong interactions of a heavy quark become independent 
of its mass and spin. 

Consider now the operators appearing at order l/mq in the effective La- 
grangian (|5l|). They are easiest to identify in the rest frame. The first operator, 

Okin = K {iD^f K ^ -^7^ K {iD f K , (56) 
2toq 2mQ 

is the gauge-covariant extension of the kinetic energy arising from the off- 
shell residual motion of the heavy quark. The second operator is the non- 
abelian analogue of the Pauli interaction, which describes the chromo-magnetic 
coupling of the heavy-quark spin to the gluon field: 

Cmag = ^M'' ^'"^ ~^ ~ ~ S- Be K ■ (57) 

4mQ rnq 



Here S is the spin operator defined in (js^), and B\ — —^e^^^G^^ are the 
components of the chromo-magnetic field. The chromo-magnetic interaction 
is a relativistic effect, which scales like l/mq. This is the origin of the heavy- 
quark spin symmetry. 



3.2 Wave-Function Renormalization of the Heavy-Quark Field in the HQET 

Besides being an effective theory for the strong interactions of heavy quarks 
with light degrees of freedom, the HQET is a consistent, renormalizable (order 
by order in l/mq) quantum field theory in its own right. In particular, it 
provides a framework for calculating radiative corrections. We shall discuss as 
an illustration the wave- function renormalization of the heavy-quark field hy. 

In quantum field theory, the parameters and fields of the Lagrangian have 
no direct physical significance. They have to be renormalized before they 
can be related to observable quantities. In an intermediate step the theory 
has to be regularized. The most convenient regularization scheme in QCD 
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is dimensional regularization in which the dimension of space-time is 

analyticaUy continued to = 4 — 2e, with e being infinitesimal. Loop integrals 
that are logarithmically divergent in four dimensions become finite for e > 0. 
From the fact that the action S — J d^x C{x) is dimensionless, one can derive 
the mass dimensions of the fields and parameters of the theory. For instance, 
one finds that the "bare" coupling constant a^^'^'^ is no longer dimensionless if 
D ^ A: dim[a^^'"^] — 2e. In a renormalizable theory, it is possible to rewrite 
the Lagrangian in terms of renormalized quantities in such a way that Green's 
functions of the renormalized fields remain finite as e — )■ 0. For QCD, one 
introduces renormalized quantities by Q^'"''' = Zg^^Q'™, A^''''' = zY^A'™, 
Q,baro _ /i^'Z^a"", etc., where /i is an arbitrary mass scale introduced to 
render the renormalized coupling constant dimensionless. Similarly, in the 
HQET one defines the renormalized heavy-quark field by h^'^'^'^ = zj^^^ h]^^. 
From now on, the superscript "ren" will be omitted. 



Figure 8: One-loop self-energy — iE(t) ■ fc) of a heavy quark in the HQET. 

In the minimal subtraction (ms) scheme, can be computed from the 
1/e pole in the heavy-quark self-energy using 

1 - = ipole of . (58) 

e ov ■ k 

As long as w • fc < 0, the self-energy is infrared finite and real. The result 
is gauge-dependent, however. Evaluating the diagram shown in Fig. ^ in the 
Feynman gauge, we obtain at one-loop order 



S(w • fc) = -ig^ tJa 



d^t 1 



{2n)D {t^ + iT^)[v ■ {t + k) + iT]] 
= -2iCFgi / dA 



d^t 1 
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where Cp — 4/3 is a colour factor, A is a dimensionful Fcynman parameter, 
and ljj = —2v ■ k > Q acts as an infrared cutoff. A straightforward calculation 
leads to 



- ^ r(i + f-^ /dz (1 - 



1 

dv ■ k TT \47r/i^ / J 



.2^r,2,r(i-,(^)-. (00, 

where we have substituted \ ^ lj {\ — z)/ z. From an expansion around e = 0, 
we obtain 

= l + (01, 

This result was first derived by Politzer and Wisexil. In the meantime, the 
calculation was also done at the two- loop order E3~E3. 

3.3 The Residual Mass Term and the Definition of the Heavy- Quark Mass 

The choice of the expansion parameter in the HQET, i.e. the definition of the 
heavy-quark mass t?t,q, deserves some comments. In the derivation presented 
earlier in this section, we chose mq to be the "mass in the Lagrangian" , and 
using this parameter in the phase redefinition in (^) we obtained the effective 
Lagrangian (|5^), in which the heavy-quark mass no longer appears. However, 
this treatment has its subtleties. The symmetries of the HQET allow a "resid- 
ual mass term" dm, for the heavy quark, provided that 5m is of order Aqcd 
and is the same for all heavy-quark flavours. Even if we arrange that such a 
term is not present at the tree level, it will in general be induced by quantum 
corrections. (This is unavoidable if the theory is regulated with a dimensionful 
cutoff.) Therefore, instead of (p2) we should write the effective Lagrangian in 



the more general formi 



3d 



K{x) = e*™«"-" P+ g(x) 
=> C,oo = hyiv-Dhy — Smhyhy . (62) 

If we redefine the expansion parameter according to mg mq + Am, the 
residual mass changes in the opposite way: 6m — > Sm — Am. This implies 
that there is a unique choice of the expansion parameter such that 5m — 0. 
Requiring Sm = 0, as it is usually done implicitly in the HQET, defineSpa,heavy- 
quark mass, which in perturbation theory coincides with the pole massE3. This, 
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in turn, defines for each heavy hadron a parameter A (sometimes called the 
"binding energy" ) through 

A {niH - toq) . (63) 

rriQ — >oo 

If one prefers to work with another choice of the expansion parameter, the 
values of non-perturbative parameters such as A change, but at the same time 
one has to include the residual mass term in the HQET Lagrangian. It can be 
shown that the various parameters that depend on the definition of mg enter 
the predictions for all physical observables in such_a way that the results are 
independent of which particular choice one adopts lj. 

There is one more subtlety hidden in the above discussion. The quantities 
mg, A and Sm are non-perturbative parameters of the HQET, which have a 
similar status as the vacuum condensates in QCD phenomenology EJ. These 
parameters cannot be defined unambiguously in perturbation theory. The rea- 
son lies in the divergent behaviour of perturbative expansions in large orders, 
which is associated with the existence oLainaularities along the real axis in the 



Borel plane, the so-called renormalons E3~0. For instance, the perturbation 
series which relates the pole mass mg of a heavy quark to its bare mass, 

= "^q" (l + ci as(TOg) -t- C2 a'^^img) + . . . + c„ a"(TOg) + • . . | , (64) 

contains numerical coefficients c„ that gmw_as n\ for large n, rendering the 
series divergent and not Borel summable EZrEa. The best one can achieve is to 
truncate the perturbation series at the minimal term, but this leads to an un- 
avoidable arbitrariness of order Amg ^ Aqcd (the size of the minimal term) . 
This observation, which at first sight seems a serious problem for QCD phe- 
nomenology, should actually not come as a surprise. We know that because of 
confinement quarks do not appear as physical states in nature. Hence, there 
is no way to define their on-shell properties such as a pole mass. In view of 
this, it is actually remarkable that QCD perturbation theory "knows" about its 
incompleteness and indicates, through the appearance of renormalon singular- 
ities, the presence of non-perturbative effects. We must first specify a scheme 
how to truncate the QCD perturbation series before non-perturbative state- 
ments such as Sm — become meaningful, and hence before non-perturbative 
parameters such as mg and A become well-defined quantities. The actual 
values of these parameters will depend on this scheme. 

We stress that the "renormalon ambiguities" are not a conceptual problem 
for the heavy-quark expansion. In fact, it can be shown quite generpJiy that 
these ambiguities cancel in all predictions for physical observables E3. The 
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way the cancellations occur is intricate, however. The generic structure of the 
heavy-quark expansion for an observable is of the form: 

observable ~ Clasimn)] ( I + + . . . ) . (65) 

V mg J 

Here C[as{mQ)] represents a perturbative coefficient function, and A is a 
dimensionful non-perturbative parameter. The truncation of the perturba- 
tion series defining the coefficient function leads to an arbitrariness of order 
Aqcd/™q, which precisely cancels against a corresponding arbitrariness of 
order Aqcd in the definition of the non-perturbative parameter A. 

The renormalon problem poses itself when one imagines to apply pertur- 
bation theory in very high orders. In practise, the perturbative coefficients are 
known to finite order in (at best to two-loop accuracy) , and to be consistent 
one should use them in connection with the pole mass (and A etc.) defined to 
the same order. 



4 Matching and Running 

In section ^ we have discussed the first two steps in the construction of the 
HQET. Integrating out the small components in the heavy-quark fields, a non- 
local effective action was derived, which was then expanded in a series of local 
operators. The effective Lagrangian derived that way correctly reproduces the 
long-distance physics of the full theory. It does not contain the short-distance 
physics correctly, however. The reason is obvious: A heavy quark participates 
in strong interactions through its coupling to gluons. These gluons can be soft 
or hard, i.e. their virtual momenta can be small, of the order of the confinement 
scale, or large, of the order of the heavy-quark mass. But hard gluons can 
resolve the spin and flavour quantum numbers of a heavy quark. Their effects 
lead to a renormalization of the coefficients of the operators in the HQET. 

Consider, as an example, matrix elements of the vector current V = q^^pQ ■ 
In QCD this current is (partially) conserved and needs no renormalization ''^'^'^ 



Its matrix elements are free of ultraviolet divergences. Still, these matrix 
elements have a logarithmic dependence on mg from the exchange of hard 
gluons with virtual momenta of the order of the heavy-quark mass. If one goes 
over to the effective theory by taking the limit mg — > oo, these logarithms 
diverge. Consequfintly, the vector current in the effective theory does require a 
renormalizationO. Its matrix elements depend on an arbitrary renormalization 
scale /i, which separates the regions of short- and long-distance physics. If /i 
is chosen such that Aqcd <SC /it <C mg, the effective coupling constant in the 
region between ji and mg is small, and perturbation theory can be used to 
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compute the short-distance corrections. These corrections have to be added 
to the matrix elements of the effective theory, which contain the long-distance 
physics below the scale /i. Schematically, then, the relation between matrix 
elements in the full and in the effective theory is 

C {rn ^) 

(y(mQ ))qcd = C'oimQ,^j.) (Vo(^))hqet H — — (Vi(^))hqet + • . • , 

(66) 

where we have indicated that matrix elements in the full theory depend on toq, 
whereas matrix elements in the effective theory are mass- independent, but do 
depend on the renormalization scale. The Wilson coefhcients Ci{mQ, ^) are 
defined by this relation. Order by order in perturbation theory, they can 
be computed from a comparison of the matrix elements in the two theories. 
Since the effective theory is constructed to reproduce correctly the low-energy 
behaviour of the full theory, this "matching" procedure is independent of any 
long-distance physics, such as infrared singularities, non-perturbative effects, 
the nature of the external states used in the matrix elements, etc. 

The calculation of the coefficient functions in perturbation theory uses the 
powerful methods of the renormalization group. It is in principle straightfor- 
ward, yet in practice rather tedious. A comprehensive discussion of most of the 
existing calculations of short-distance corrections in the HQET can be found 
in Rcf. 8. Here, we shall discuss as an illustration the renormalization of the 
1 /mg-suppressed operators in the effective Lagrangian (|l]) . At the tree level, 
there appear two operators at order 1 /toq , which have been given in (j5^) and 
(p7|). Beyond the tree level, the coefficients of these operators may be modi- 
fied, and other operators not present at the classical level may be induced. In 
general, we thus expect 

Ci/'m = Ci,in{p) Ckin(M) + C,^a^{p) Oi„ag(M) + ^cw Operators , (67) 

where /i is the renormalization scale. But how do we calculate the Wilson 
coefficient functions, and what are the possible new operators? To extend the 



classical construction of Sec. 3.1 to include quantum corrections would be cum- 
bersome. Fortunately, there is a systematic procedure which allows us to derive 
the result in the presence of quantum effects in a rather simple and straight- 
forward way. It consists of three steps: construction of the operator basis, 
calculation of the "matching conditions" aX ji — niQ, and renormalization- 
group improvement ("running"). Below, we shall first explain these steps in 
general and then illustrate them with the particular example of Ci/^. 
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Construction of the Operator Basis 

Similar to the fields and coupling constants, in a quantum field theory any com- 
posite operator built from quark and giuon fields may require a renormalization 
beyond that of its component fields. Such operators can be divided into three 
classes: gauge-invariant operators that do not vanish by the equations of mo- 
tion (class-I) , gauge- invariant operators that vanish by the equations of motion 
(class-II), and operators which are not gauge-invariant (class-Ill). In general, 
operators with the same dimension and quantum numbers mix under renor- 
malization — Hossjever, things simplify if one works with the background field 
technique ll2iH12J, which is an elegant method for quantizing gauge theories, 
preserving explicit gauge invariance. This offers the advantage that a class- 
I operator cannot mix with classJII operators, so that only gauge-invariant 
operators need to be consideredll23. Furthermore, class-II operators are irrele- 
vant since their matrix elements vanish by the equations of motion. It it thus 
sufficient to consider class-I operators only. 

Thus, we must find a complete set of class-I operators of the right dimen- 
sion, carrying the quantum numbers allowed by the symmetries of the problem. 
In the case at hand, we are dealing with operators appearing at order l/mg 
in a strong-interaction Lagrangian, and we thus have to find dimension-five 
operators containing two heavy-quark fields of the same velocity. Moreover, 
these operators must transform as scalars under the Lorentz group. The most 
general form of such operators is 

JiyT^^iD'^'iD'' hy; T^t/ e jspi^^ I'/^Wi', Ti'^^a" I^a" } ■ (6^) 

Note that the velocity is not a dynamical quantity in the HQET and thus can 
be used to construct the basis operators. Using that hyj^ = hy hy, we 
find that there are only three possible operators: 

hyilDfhy, hy{ivDfhy, ^ hy g sC^" hy . (69) 

Since the equation of motion of the HQET is iv-D hy = 0, it follows that there 
are two class-I and one class-II operators, which we choose in the form: 

class-I: hy{iD^yhy, \hy af^^gsG^'^hy , 
class-II: hy{ivDfhy. (70) 

Besides a class-II operator, which has vanishing matrix elements between phys- 
ical states, the kinetic and chromo-magnetic operators already present at the 
tree level are thus the only operators which can appear in Li/m, even in the 
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presence of quantum corrections. Once we have found a complete basis of 
class-I operators, our next goal is to calculate their coefficient functions in 
perturbation theory. 



4-2 Matching Conditions at fi = mq 

The Wilson coefficient functions Ckin(M) and Ci„ag(A') in ( |67| ) can be obtained 
from the comparison ( "matching" ) of Green's functions in QCD with those in 
the effective theory. It is crucial that, by construction, the Wilson coefficients 
receive only short-distance contributions (see Fig. ^) and are thus insensitive 
to the properties of the external states. This ensures that once the coefficients 
have been determined by requiring that some particular Green's function(s) 
be the same in the two theories, all other Green's functions will be the same. 
Moreover, since the Wilson coefficients are infrared insensitive they are calcu- 
lable in perturbation theory, and we can perform their calculation using quark 
and gluon states rather than physical hadron states. 

In the example at hand, the coefficients Ckin(M) and Cmag(/i) can be ob- 
tained from a calculation of the Green's function of two heavy quarks and a 
background gluon field, to one-loop order in the full and in the effective the- 
ory. The relevant vertex diagrams in QCD are shown in Fig. |^. They have to 
be supplemented by the wave-function renormalization of the external quark 
lines. The background field is not renormalized. The momentum assignments 
are such that p is the outgoing momentum of the background field, and k and 
(fc — p) are the residual momenta of the heavy quarks. To order 1/mq, it is 
sufficient to keep terms linear in k or p. The quarks can be taken on shell, in 
which case v ■ k — v ■ p — 0. A subtlety which has to be taken into account is 
that the QCD spinor uq{Pq,s) is related to the spinor Uh{v, s) of the effective 
theory by 

uq{Pq,s)=(l + -^ + ...)uh{v,s), (71) 



2mQ 

where Pq — mqv + k. In the matching calculation one has to use the same 
spinors in both theories. We thus define a vertex function F'' by writing the 
amplitude as i^'^ gsAfj^ a^ji) UhX^taUh, so that at the tree level in QCD 

- (^ ^ ^^1^^ 



= .^+(^^-^)%l2^ + .... (72) 
2mq Amq 

Here the ellipses represent terms of higher order in k or p, and we have used 
that between the heavy-quark spinors 7'^ can be replaced by . 
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niQii+k mqU+k-p 

Figure 9: Diagrams for the calculation of the heavy quark-gluon vertex function in 
QCD. The background field is denoted by A. 



The contributions to the vertex function arising at the one-loop level are 
also shown in Fig. ^. They contain both abelian and non-abelian vertices. 
Since the matching calculation is insensitive to any long-distance properties 
such as the nature of the infrared regulator, it is legitimate to work with 
anyfixLfcared regularization scheme that is convenient. Following Eichten and 



nyjxLtca 

HillLa'1123, we choose to regulate both ultraviolet and infrared divergences using 
dimensional regularization. Moreover, we expand the resulting expressions for 
the Feynman amplitudes to linear order in the external momenta and then set 
the external momenta to zero inside the loop integrals. Then the only mass 
scale remaining is the heavy-quark mass. In the MS scheme, the result for the 
one- loop contribution to the QCD vertex function islB 

r(^CB.-^^f-^.ln^ + C. + C.V (73) 



where Cp = ^{N^ ~ l)/-^c = 4/3 and Ca = A^c = 3 are the eigenvalues of the 
quadratic Casimir operator in the fundamental and the adjoint representations. 

Now comes a clue: if dimensional regularization is used to regulate both ul- 
traviolet and infrared singularities, all loop integrals in the HQET are no-scale 
integrals (after a power of the external momenta has been factored out) and 
vanish! So only the tree-level matrix elements of the HQET operators in the 
effective Lagrangian multiplied by their Wilson coefficient functions remain. 
This is why dimensional regularization is superb for matching calculations. 
The result is 

rJ^QET =^' + ^^^^ + Cmag(M) + ■■■■ (74) 

Requiring that the vertex functions be the same in the full and in the effective 
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theory, we find (in the MS scheme) 



Ckin(A^) = l, a,ag(M)-l + 77^ -CAln^ + C^ + C;^ . (75) 



The fact that the kinetic operator is not renormalized is not an accident, but 
follows from an invariance of the HQET under small redefinitions of the velocity 
used in the construction of the effective Lagrangian. Clearly, the predictions 
of the HQET should not depend on whether v is taken to be the velocity 
of the hadron containing the heavy quark, the velocity of the heavy quark 
itself, or some other velocity differing from the hadron velocity by an amount 
of order Aqcd /wq • This so-called reparametrization invaxiaxLce implies that 
Ckin(/^) = 1 must hold to all orders in perturbation theory E£3'E£3. 

In the next paragraph, we will see that the scale dependence predicted by 
the one-loop result quoted above cannot be trusted if ^ ™q; however, what 
can be obtained from the matching calculation are the values of the coefficient 
functions at the matching scale fi — toq as well as their logarithmic derivatives. 
For the coefficient of the chromo-magnetic operator, we find: 



4-. 3 Renormalization- Group Evolution 

The one-loop calculation presented above allows us to derive expressions for 
the Wilson coefficient functions provided that mg/fi = 0(1). In practical 
applications of effective field theories, one is however often interested in the 
case where there is a large ratio of mass scales. After all, an effective theory 
is constructed to separate the physics on two very different energy scales. In 
such a situation, the coefficient functions contain large logarithms of the type 
[cks ln(mQ//x)]", which must be summed to all orders in perturbation theory. 
This is achieved by using the powerful machinery of the renormalization group. 

For a set {Oi} of n class-I operators that mix under renormalization, one 
defines an n x n matrix of renormalization factors Zij by O^^'''^ = Zij Oj (fj.) , 
such that the matrix elements of the renormalized operators Oj{fj,) remain 
finite as e — > 0. In contrast to the bare operators, the renormalized ones 
depend on the subtraction scale via the /i dependence of Zij : 



Cinag("lQ) 
dlnCmag(M) 

dln/i 



1 + (Ca + Cf) 



2tt 




(76) 



ibaic 



(77) 
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where 



d/i 



(78) 



are called the anomalous dimensions. That under a change of the renormal- 
ization scale the operators mix among themselves follows from the fact that 
the basis of operators is complete. It is convenient to introduce a compact 
matrix notation, in which 0(/i) is the vector of renormalized operators, Z is 
the matrix of renormalization factors, and 7 denotes the anomalous dimension 
matrix. Then the scale dependence of the renormalized operators is controlled 
by the renormalization-group equation (RGE) 



(79) 



/^^+7)O(a*)=0. 



In the MS scheme, the matrix Z obeys an expansion of the form 



Z = 1 



oc ^ 
k=l 



Zh (a.. 



(80) 
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and by requiring that the anomalous dimensions in (|7^) be finite as e ~^J^ 
finds that 7 can be computed in terms of the coefficient of the 1/e pole 

dZiias) 



one 



7 



-2a., 



da. 



(81) 



The same relation holds in the MS scheme. 

From (|79| ) and the fact that the product Ci (fi) Oi (/i) must be /i indepen- 
dent, we derive the RGE satisfied by the coefhcient functions. It reads 



d/i 



(82) 



where we have collected the coefficients into a vector C{ii). In general, the 
Wilson coefficients can depend on both explicitly or implicitly through the 
running coupling. We thus have 



_d_ _ _9_ 
d/i ^ 



d 



das{p) ' 



where the (3 function 



in 



as 
An 



(83) 



(84) 
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e-one- 



iic 



describes the scale dependence of the renormahzed couphng constant— II 
and two- loop coefficients are scheme independent and arc given bycSiij 

Po = Y 3 ' 

Pi = -Cl- i-CA + ^CFjTpnf, (85) 

where n/ is the number of light quark flavours, and Tp — 1/2 is the normal- 
ization of the SU(3) generators in the fundamental representation: tr(taib) = 
Tp Sab- It is now straightforward to obtain a formal solution of the RGE. It 
reads 

Cifi)^U{ti,mQ)C{mQ), (86) 
with the evolution matrix EiUlii 

C/(M,mQ) =r„ exp J da^^. (87) 

Here "Tq," means an ordering in the coupling constant such that the couplings 
increase from right to left (for /i < toq). This is necessary since, in general, 
the anomalous dimension matrices at different values of a, do not commute. 



Eq. IpA) can be solved perturbatively by expanding the f3 function (84) and the 



anomalous dimension matrix in powers of the renormalized coupling constant: 

^ . as f as 

Here we shall only discuss the important case of a single coefficient function, 
or equivalently, when there is no operator mixing. Then the matrix 7 reduces 
to a number, and the evolution is described by a function U (/i, toq), for which 
the perturbative solution of (|8^) at next-to-leading order yields 

with 

70 c, 71 loPi /nn^ 

" = ' = Wo-M- ^'"^ 

The theoretical framework discussed here is called "renormalization-group 
(RG) improved perturbation theory" . In the expression for the evolution func- 
tion J7(^, mg), there are no large logarithms of the form ashiirnq / fj,) left. 



5 + , (89) 
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They are all contained in the ratio of the running couplings evaluated at the 
scales niQ and /i. Thus, RG-improved perturbation theory provides the optimal 
method to bridge wide energy intervals. (As a side remark, we note that the 
same technique is used to control the evolution of gauge couplings and running 
mass parameters from low energies up to very high energy scales characteristic 
of grand unified theories.) The terms shown explicitly in (|89| ) correspond to the 
so-called next-to-leading order (NLO) in RG-improved perturbation theory. In 
this approximation, the leading and subleading logarithms [as \n(mQ / (!)]"■ and 
as[Q!s ln(TOQ/^)]" are summed correctly to all orders in perturbation theory. 
To achieve this, it is necessary to calculate the two-loop coefficient 71 of the 
anomalous dimension. When 71 is not known, it is only possible to evaluate 
the evolution function in the so-called leading logarithmic order (LO), in which 

UM,,mc,)^(^y . (91) 



as (a*) 

This still sums the leading logarithms to all orders, but does not contain the 
non-logarithmic terms of order as- 

To complete the calculation of the RG-improved coefficient function, the 
evolution function U{fi,mQ) must be combined with the initial condition for 
the Wilson coefficient at the high energy scale /i = toq, as shown in (^6|). 
If the operator under consideration is present at the tree level, the matching 
condition can be written in the form 

C(mQ) = l + ci^^ + ..., (92) 

where ci is obtained from a one-loop calculation. To obtain a consistent (i.e. 
renormalization-scheme independent) result at next-to-leading order, we have 
to combine the one-loop matching condition with the expression for J7(/i,TOQ) 
given in (|8^). This requires the calculation of the two- loop anomalous dimen- 
sion. The result is 

In this expression, the terms involviiig—the coupling constant asimq) are 



renormalization-scheme independent The exponent a involves only the 



one-loop coefficients 70 and (3o and is scheme independent by itself. For the 
coefficient {S + ci) of the next-to- leading term things are more complicated, 
however. The one- loop matching coefficient ci, the two- loop anomalous di- 
mension 71, and the QCD scale parameter Aqcd in the expression for the 
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running coupling constant arc all scheme dependent, but they conspire to give 
asimg) a scheme-independent coefficient. On the other hand, the coefficient S 
ofasi/J.) does depend on the renormalization procedure. This is not a surprise; 
only when the /Lt-dependent terms in the Wilson coefficient are combined with 
the /x-dependent matrix elements of the renormalized operator one can expect 
to obtain a scheme-independent result. For this reason, it is sometimes useful 
to factorize the solution (|9^) in the form C(/i) = C{mQ)K{fi), where C{mQ) 
is RG-invariant and contains all dependence on the large mass scale toq. The 
scheme-dependent function K{^) can be used to define a RG-invariant renor- 
malized operator: O = K {/j,) O (fj,) . At next-to-leading order, we obtain: 

diruQ) = [asimQr |l + ^^(5 + Ci)} , 

a=K(A.r''|l-^5|o(A*). (94) 

Let us finally apply this formalism to the operators appearing at order 
1/toq in the effective Lagrangian of the HQET. The fact that reparametriza- 
tion invariance ensures that the kinetic operator is not renormalized implies 
that the 2x2 anomalous dimension matrix for the operators Okin and Omag 
is diagonal and of the form 

7 = -f^^s ) ■ (95) 

The one-loop coefficient of the anomalous dimension of the chromo-magnetic 
operator, together with the one-loop matching coefficient, can be obtained 
from (|7|): 

70"^^^^ = 2Ca , = 2{Ca + Cf) . (96) 

The result for 7™^^ can also be obtained in a simpler way bi computing only 
the 1/e poles in the matrix elements of the bare operators E3. Unfortunately, 
the two-loop coefficient 7™*^^ is not yet know n .f^ This means that the coefficient 
Smag in the next-to-leading order solution (^^ is still unknown. 

4-4 Renormalization of Heavy- Quark Currents 

As a final example, we discuss the renormalization of local current operators 
involving two heavy-quark fields. This case is of particular importance for 

''This is one of the few cases of interest where an anomalous dimension in the HQET is not 
yet known to two-loop order. If you feel strong enough, you are invited to try the calculation 
of 7^''! 
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phenomenology, as the weak current for b c£i> transitions is of this form. 
In the HQET, the relevant current operators contain two heavy-quark fields at 
different velocity and are thus of the form hyiT hy (it does not matter whether 
the two fields have the same flavour), where F is some Dirac matrix, whose 



structure is irrelevant to our discussion. We have discussed in Sec. 2.3 that the 
matrix elements of such operators between meson states are proportional to 
the universal Isgur-Wise form factor ^(v ■ v'). We shall now derive, in leading 
logarithmic order, the Wilson coefhcient function that relates the QCD current 
operators with their HQET counterparts rcnormalizcd at the scale fj, toq. 




Figure 10: One-loop vertex diagram arising in the calculation of the anomalous 
dimension of heavy-quark currents. The external current changes the heavy-quark 
velocity from v to v' . 

To this end, we need to calculate, using dimensional regularization, the 
1/e pole in the matrix element of the bare current operator between quark 
states. The relevant vertex diagram is shown in Fig. Since in the effective 
theory the coupling of a heavy quark to a gluon does not involve a 7 matrix, 
it is easy to see that to all orders in perturbation theory the operator hy'T hy 
is renormalized multiplicatively and irrespective of its Dirac structure. The 
extraction of the one-loop ultraviolet divergence can be done in a few lines. In 
the Feynman gauge, the value of the vertex diagram is (omitting the quark 
spinors) : 



'Aig^tJaV ■ v'T 



d^i 1 



(27r)^ (t2 -^iri)(^v .t + iri)[v' ■ t + irj) 



d^t 1 







(27r)^ [t2 + 2{pv + Xv') - t + iri] 



^The final result for the anomalous dimension is gauge independent. 
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oc oo 



= _^r(l + e)..t;'r(4Vr fdX /dp-— --J-^-^. 

TT J J (p^ + x^ + 2v-v'pxy^ 



(97) 

Defining a new variable z — p/X, and introducing an arbitrary infrared cutoff 
S, we can rewrite tlie double integral in the form: 

oo oo 

dX 7— / dz z-r- = r(w) + finite terms, (98) 

(A2 + 52)i+«^7 (l + z2 + 2wz)^+^ 2e ^ ^ ' ^ ^ 

u 

where w = v ■ v' , and 

oo 

r{w) = /dz ^— — = , I In fw + ^w"^ - l] . (99) 

J l + z^ + 2wz - 1 V y ^ ^ 



Hence, the 1/e pole of the vertex diagram is given by 

Cpa. 



27re 



ru;r(u;). (100) 



To obtain the renormalization constant Z^h of the bare current operator, we 
have to add a contribution 

from the wave-function renormalization of the heavy-quark fields, where 
has been given in (|6l|). The result is 

C a 

Zhh = 1 ^ — - [wr{w) — 1] + finite terms. (102) 

ZTTC 

By means of the relation (^l|), we derive from this the one- loop coefficient of 
the anomalous dimension of heavy-quark currents in the HQET. This is_.the 
famous velocity-dependent anomalous dimension obtained by Falk et al. l3: 

-fl;^{w) = ACf [w r{w) - 1] , (103) 

In the zero-recoil limit, i.e. for w — 1, the heavy-quark currents are the sym- 
metry currents of the spin-flavour symmetry, and as such they are not renor- 
malized, since the associated charges are conserved. This implies that 

7'''^(1) = (104) 
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to all orders in perturbation theory. This constraint is satisfied by the one-loop 
result in (103), since r(l) — 1. 

Since in the effective theory the velocity of a heavy quark is conserved iy 
the strong interactions, the heavy quark can be described by a Wilson linetZl. 
An external current can instantaneously change the velocity, resulting in a 
kink of that line. It is well-known that such cusps lead to singular behaviour. 
The renormalization of cusp singularities_Df Wilson lines was investigated in 



detail by Korchemsky and Radyushkin Eij already in 1987, prior to the de- 
velopment of the HQET. In particular, they calculated the one- and two-loop 
coefficients of the so-called cusp anomalous dimension ^'^^^^{ip) as a function of 
the hyperbolic cusp angle ip. But this anomalous dimension is precisely that of 
heavy-quark currents with the identification cosh(p = w. Later, the result 
for 75'^(u'), wiiicjh we will not present here, has been confirmed in the context 
of the HQETL 



lie 




Figure 11: Velocity dependence of the Wilson coefficient ChhCwiA*): evaluated for 
as{ii) I asiraq) = 2 and Uf = 3. 

In leading logarithmic order, the expansion of heavy-quark currents in the 
HQET takes the form 



QTQ ^ Ci,h{w,iJ.) hyiThv + 0{l/mQ) , 



(105) 



where 

Chh(w, m) 



2Cf 



[wr{w) - 1] . (106) 



The velocity dependence of the Wilson coefficient is illustrated in Fig. |ri|. 
Physical decay or scattering amplitudes are proportional to the product 



Chh(w,/i)^(w,/i) , 



(107) 
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where £,{w, /i) is the renormaUzed Isgur-Wise function, which satisfies ^(1, /i) = 
1. The fact that ahh(l) = imphes that RG effects respect the normalization 
of form factors at zero recoil. 

We like to finish this discussion with a curious remark, which illustrates 
the power of RG methods. For values w = 0(1), the leading-order expres- 
sion for the coefficient Chh('ii',A') in ( |lO(j ) contains all large logarithms of the 
type [tts ln(mQ//L<)]". For very large values of w, however, we have wr{w) — > 
ln(2w) — 1, so that in the RG improvement of Ch\i{w,^) we have resummed 
terms of the form [as In w ln(mQ //i)]". These are the well-known Sudakov dou- 
ble logarithms, which arise from the emission of gluon bremsstrahlungs during 
the scattering of the heavy quarks. This effect leads to a fractional power- like 
damping of the transition form factors at large recoil, which-adds to the "soft" 
suppression contained in the Isgur-Wise form factor itself till. Explicitly, we 
obtain 

a.K.)->(f)\ .= ^ln^. (108) 
\2w/ Po as[mQ) 



5 Concluding Remarks 

We have presented an introduction to heavy-quark symmetry, the heavy-quark 
effective theory and the 1/toq expansion, which provide the modern theoret- 
ical tools to perform quantitative calculations in heavy-flavour physics. Our 
hope was to convince the reader that heavy-flavour physics is a rich and di- 
verse area of research, which is at present characterized by a fruitful interplay 
between theory and experiments. This has led to many significant discoveries 
and developments on both sides. Heavy-quark physics has the potential to 
determine many important parameters of the electroweak theory and to test 
the Standard Model at low energies. At the same time, it provides an ideal 
laboratory to study the nature of non-perturbative phenomena in QCD, still 
one of the least understood properties of the Standard Model. 

Let us finish with a somewhat philosophical remark: At this school, we 
have heard a lot about exciting new developments related to dualities, which 
relate apparently very different theories to each other. So are electric, weak- 
coupling phenomena in one theory dual to magnetic, strong-coupling phenom- 
ena in another theory. Some people argue quite convincingly that duality seems 
to be everywhere in nature, and consequently there are no really difficult ques- 
tions in physics; very difficult problems become trivial when approached from a 
different, dual point of view. There are, however, "moderately difficult" prob- 
lems in physics, which are "self-dual" . It is the author's opinion that real- world 
(i.e. non-supersymmetric) QCD at hadronic energies belongs to this category. 
Having said this, we conclude that heavy-quark effective theory provides a 
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powerful tool to tackle the "moderately difficult" problems of heavy-flavour 
physics. 
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